Low-energy limit of the radiative dipole strength in nuclei 
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We explain the low-energy anomaly reported in several experimental studies of the radiative dipole 
strength functions in medium-mass nuclei. These strength functions at very low gamma- energies 
correspond to the gamma-transitions between very close nuclear excited states in the quasicontin- 
uum. In terms of the thermal mean-field, the low-energy enhancement of the strength functions in 
highly-excited compound nuclei is explained by nucleonic transitions from the thermally unblocked 
single-quasiparticle states to the single-(quasi)particle continuum. This result is obtained within 
the finite-temperature quasiparticle random phase approximation in the coordinate space with ex- 
act treatment of the single-particle continuum and exactly eliminated spurious translational mode. 
The case of radiative dipole strength functions at the nuclear excitation energies typical for the 
thermal neutron capture is illustrated for ^*Mo and ^'*''Nd in comparison to available data. 

PACS numbers: 21.10.Pc, 21.60.Jz, 25.40.Lw, 27.60.+j 
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Experimental and theoretical studies of the nuclear 
low-energy electric dipole response remain among the 
challenges of the modern nuclear structure physics and 
attract an increasing interest because of its astrophysical 
impact. Radiative strength (7-strength) at low energies 
may enhance the neutron capture rates in the r-process of 
nucleosynthesis [l|, ^ with a considerable influence on el- 
emental abundance distributions. One of the key phases 
of the r-process nucleosynthesis is capture of a thermal 
neutron with the subsequent 7-decay of the compound 
nucleus. The typical neutron energy in the astrophysi- 
cal plasma is about 100 keV. Therefore, the description 
of 7-emission spectra of a compound nucleus with ex- 
citation energies of the order of the neutron separation 
energy is the central problem. Hauser-Feshbach model is 
a standard tool for calculations of the radiative neutron 
capture cross sections [3j. Formally, this model includes 
all possible decay channels via transmission coefficients. 
In the gamma-decay channel the corresponding coeffi- 
cient is determined by the radiative strength function 
which is usually calculated by one of the phenomeno- 
logical parameterizations • However, in more recent 
works [1, 2, 7] it has been shown that for the most impor- 
tant electric dipole strength these simple models are not 
sufficient because they do not account for structural de- 
tails of the strength at the neutron threshold. Sensitivity 
of the stellar reaction rates to these details emphasizes 
the importance of their studies within microscopic self- 
consistent models. 

Another key ingredient for the Hauser-Feshbach cal- 
culations is the Brink-Axel hypothesis [§] stating that 
the 7-strength does not depend on the nuclear excitation 
energy, in particular, it is the same for excited and non- 
excited nuclei. Supposedly true for the giant resonances 
and for the soft modes like pygmy dipole resonance, this 
hypothesis is, however, violated for the lowest transition 
energies. For instance, non-zero strength is systemati- 



cally observed at very low gamma-energies . Radiative 
strength functions extracted from various measurements 
[l0l - [l4| show an upbend at Ej < 3 MeV in light nuclei 
of Fe-Mo mass region. Studies of Ref. have revealed 
that this phenomenon, occuring in various astrophysical 
sites, can have a significant impact on their elemental 
abundances. Phenomenological approaches approximate 
the low-energy 7-strength by the tail of the giant dipole 
resonance with a temperature-dependent width. This is, 
however, not justified, because the low-energy 7-strength 
originates from underlying physics which is completely 
different from the giant vibrational motion. Modern 
microscopic theories have excellent tools for computing 
probabilities of transitions between the nuclear ground 
state and excited states, but have common problems to 
describe transitions between excited states. The gen- 
eral many-body techniques like Green function formalism 
p^.[T7j. can be applied to 7-emission and 7- absorption in 
excited states of compound type if it is approximated by 
a semi-statistical model like a finite-temperature mean- 
field. In such a case, the transitions are described by 
the finite-temperature version of the random-phase ap- 
proximation and its extensions. There exist formulations 
within discrete model spaces [T8l-l20| and models with ex- 
act treatment of the single-particle continuum [2l| - [2^ . 

In this Letter, we explain the mechanism for the en- 
hancement of the low-frequency dipole 7-transitions be- 
tween the nuclear excited states in the quasicontinuum. 
We show, for the first time, that this phenomenon can be 
quantitatively described in terms of a microscopic many- 
body approach built on the thermal mean-field descrip- 
tion of the compound nucleus. Exact treatment of single- 
particle continuum at finite temperature and exact elimi- 
nation of the center-of-mass motion are the two essential 
ingredients for understanding these dipole 7-transitions 
with frequencies < 3-4 MeV. First of all, we show 
that the thermal mean field model provides a reasonable 
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approach to the compound nucleus. At the same time, 
this model is simple enough to allow a straightforward 
generalization of the many-body methods for the case 
of a compound nucleus, in terms of finite temperature 
corresponding to the nuclear excitation energy. The ap- 
plicability of the finite-temperature mean-field approach 
to compound nuclei is usually justified by the fact that 
the product of the canonical distribution and the nuclear 
level density entering the integrals for thermal averages 
is a sharp peak approximating the delta-function of the 
micro-canonical distribution. Indeed, the nuclear level 
density has an exponential dependence on energy which, 
together with the canonical distribution, gives a bell- 
shaped function of energy. But for real medium-mass 
and heavy nuclei the typical width of this function is 
estimated by several MeV, within the conditions corre- 
sponding to the thermal neutron capture. However, as 
it is shown below, the dependence of the nuclear excita- 
tion energy on the temperature parameter in the thermal 
mean-field approach is consistent with their phenomeno- 
logical dependence and with empirical level densities [1^ . 
This fact is used as an alternative justification for em- 
ploying the finite-temperature mean-field approach for 
compound nuclei. 

The general concept of the finite-temperature mean- 
field theory [ll[ll,[2| is based on the variational princi- 
ple of maximum entropy minimizing the thermodynami- 
cal potential 
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n{X,T) = E - XN ~TS, 



(1) 



with the Lagrange multipliers A and T determined by the 
average energy E, particle number N and the entropy S. 
These quantities are the thermal averages involving the 
generalized one-body density operator TZ: 

5 = -fcTr(7^1r^7^), E^TviWH), Tr(7^AA), (2) 

where Ti is the nuclear Hamiltonian, M is the particle 
number operator, and k is Boltzmann constant. Varying 
the Eq. ([l}, one can determine the density operator TZ 
with the unity trace: 
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The form of the nuclear energy density functional E[TZ] 
is a subiect of extensive studies, see, for instance, a re- 
view [27[. Density functional theory is supposed, in prin- 
ciple, to provide an exact description of many-body prob- 
lem [2^], if the exact density functional is known, how- 
ever, nuclear physics is still far from a microscopic deriva- 
tion of such a functional. A covariant density functional 
theory (CDFT) [2^ is one of the most successful exam- 
ples since it combines a density dependence through a 
non-linear coupling between the meson fields [30| with 
pairing correlations based on an effective interaction of 
finite range. The Lorentz invariance not only allows one 
to describe the nuclear spin-orbit coupling, which has an 
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FIG. 1: The dependence of the excitation energy of the ther- 
mally excited nuclei on the temperature parameter within the 
TRMF model (solid curves), in comparison with the parabolic 
fit (dashed curves). 

essential influence on the underlying shell structure, in a 
consistent way, but it also puts stringent restrictions on 
the number of parameters in the corresponding function- 
al without reducing the quality of the agreement with 
experimental data. Although the CDFT cannot provide 
an exact treatment of the full nuclear dynamics, it is a 
very good starting point for implementation of correla- 
tions which are taken in the mean-field only in a static ap- 
proximation. The CDFT is very convenient for our pur- 
pose providing the self-consistent mean-field description 
where the nucleonic potentials are calculated from the nu- 
cleonic densities. Fig. [T] displays the dependence of the 
excitation energy E* for selected nuclei on the tempera- 
ture parameter T obtained within the thermal relativis- 
tic mean-field (TRMF) model with NL3 parametrization 
jsij l . The dependence E* (T) shows features of the back- 
shifted Fermi gas (FG) distribution E*{T) = 5 + aT'^, 
therefore we have found reasonable to make the parabolic 
fits. The values of the level density parameter auMF de- 
termined from the TRMF model are rather close to the 
empirical values apc ~ A/8 MeV~^. Thus, the temper- 
ature of an excited compound nucleus can be obtained 
from its excitation energy by the self-consistent finite- 
temperature mean-field calculations. In this work, we 
consider excitation energies E* around the neutron sep- 
aration energies Sn of stable medium-mass and heavy 
nuclei that corresponds, approximately, to the tempera- 
tures around and below 1.5 MeV. At such temperatures 
the mean-field single-particle energies and wave functions 
remain almost unchanged. 

7-emission and 7-absorption of the compound nuclear 
system are described as an interaction of the nucleus with 
a sufficiently weak external electromagnetic field P oscil- 
lating with some frequency uj. The interaction with such 
a field causes small amplitude nuclear oscillations around 
the static equilibrium, so that the total density matrix TZ 
has an oscillating term, in addition to the static thermal 
mean-field part 7^": 



TZ{t) ^TZ°+ [STZe-"^' + h.c] . 



(4) 



0123456789 



FIG. 2: Schematic picture of the possible lowest-energy single- 
quasiparticle transitions from the thermally unblocked states 
in an excited compound (left) and from the "frozen" ones in 
the ground-state nucleus (right). 

Variation STZ of the density matrix TZ in the external 
field P obeys, in the local approximation, the following 
integral equation: 



dx'dx"Aix, x';uj, T)F{x', x")Snix";uj, T), (5) 



where a; is a multi-index x ~ {r, s, r, x} of spatial coor- 
dinate r, spin s, isospin r and component in the quasi- 
particle space x- F{x,x') is the effective nucleon-nucleon 
interaction, A{x,x';u!,T) is the two-quasiparticle propa- 
gator in the nuclear medium at finite-temperature and 
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j dx'A{x,x']w,T)P{x'). 
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The propagator A{x,x' ]uj,T) is the key quantity and 
ideally has to include all the in-medium and surface ef- 
fects. In the first approximation we calculate it within 
the thermal continuum quasiparticle random phase ap- 
proximation (TCQRPA) in terms of the Matsubara tem- 
perature Green functions [l^, The full expression 
of the TCQRPA propagator in the coordinate space is 
presented in [2^. As the effects of finite temperature 
and continuum are supposed to give the leading order 
contribution to the 7-strength at \uj\ — i> 0, at this step 
we use a non-relativistic approach similar to that of Ref. 
[23} . Thus, in our first illustrative calculations the nu- 
cleonic wave functions result from the phenomenological 
mean field and the effective nucleon-nucleon interaction 
F{x,x') has the Landau-Migdal ansatz. The propaga- 
tor consists of the discrete and continuum parts. The 
discrete part describes transitions between the single- 
quasiparticle states in the discrete spectrum, and the 
continuum part describes transitions from the discrete 
spectrum states to the continuum. Dashed and solid 
arrows in Fig. [2] show the low-frequency transitions of 
both kinds, respectively. The effective occupation prob- 
ability distribution hi{Ei,T) has much larger diffuseness 
at T = Ts = T{Sn) than at T=0, being the following 
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FIG. 3: The El 7-strength function for the thermally excited 
state of ^*Mo at the temperature parameter T = Ts = T{Sn) 
It is obtained within the TCQRPA, compared to data [ij] and 
to the 7-strength for the ground state (T=0, dashed curve). 
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FIG. 4: Same as in Fig. O but for ^^^Nd, compared to data 
from Ref. [H. 

product: h.,{Ei,T) = v'i[T){l ~ n.i{Ei,T)), where are 
the occupation numbers of Bogoliubov quasiparticles. 
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exp(S,(T)/fcr)' 
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and Ei are the eigenvalues of the single-particle Hamil- 
tonian. The radiative dipole strength function is deter- 
mined by the quantity STZ through its convolution with 
the electromagnetic dipole operator Pei- 



fEi{Ey,T) — - 



21{hc) 



w = i A, A ^- 0. Formally, Eq. 
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corresponds to 
7-absorption, and 7-emission strength can be calculated 
for the "final temperature" T/ = y/{E* - 5- E^)/a [12]. 
However, for E-^ < 3-4 MeV the 7-absorption and 7- 
emission strength functions almost coincide. Their dif- 
ferences will be discussed elsewhere. Figs. [3] and |3] dis- 
play the radiative dipole strength functions in ^''Mo and 
^'^''Nd calculated within the TCQRPA at finite and zero 
temperatures. The values of the finite temperature pa- 
rameter correspond, in each case, to the S„ value and. 



4 



hence, to the typical compound nucleus formed after the 
thermal neutron capture. The finite imaginary part A 
of the energy variable is taken equal to 200 keV and 150 
keV, respectively, in accordance to the energy resolution 
of the experimental data, in order to make a correct com- 
parison. While the dashed curves (T=0, ground state) 
show at low 7-energies the tails of the higher-energy tran- 
sitions between discrete states (blue dashed arrows in Fig. 
[2]) only due to the non-zero value of A, the solid curves 
give in addition the pure thermal continuum strength 
which remains finite at A = 0. The Figs. [3] and H] il- 
lustrate the behavior of the dipole 7-strength which is 
typical for nuclei in these two mass regions: the upbend 
of the 7-strength appears as a typical feature of the re- 
sponse of medium-mass nuclei in Fe-Mo region while in 
heavier nuclei like ^^^Nd the 7-strength comes out rather 
flat at Ej ^ 0. 

Fig- E] gives a qualitative interpretation for the low- 
energy enhancement of the 7-strength. It is clearly seen 
that transitions from the thermally unblocked states of 
the single-particle spectrum to the continuum (solid ar- 
row in Fig. m left part T>0) form solely the 7-strength 
function at very low transition frequency E^. Such type 
of transitions is not possible in the ground state (T=0) 
where the lowest-energy transitions to the continuum 
have much higher energies (solid arrow in Fig. [51 right 
part). This schematic picture also explains why the up- 
bend of the 7-strength becomes less pronounced in heav- 
ier systems: the temperature parameter corresponding to 



the neutron threshold Sn is smaller since the threshold 
itself is typically lower than in lighter systems. Accord- 
ingly, the effective thermal occupation probabilities n{Ei) 
for the single-nucleon states closest to the continuum de- 
crease that, in turn, decreases the transition probabili- 
ties. 

Summarizing, we give a theoretical interpretation of 
the low-energy anomaly in the behavior of the radiative 
dipole strength in medium-mass and heavy nuclei. We 
have shown that a microscopic approach to nuclear re- 
sponse with coupling to the continuum and exactly elim- 
inated center-of-mass motion, based on the statistical de- 
scription of the compound nucleus, gives a very good ap- 
proximation to the low-energy 7-strength already on the 
level of the two-quasiparticle configurations. Application 
to electric dipole response explains the systematic low- 
energy enhancement of the 7-strength on the microscopic 
level. Thus, it is shown that microscopic nuclear many- 
body theory can be brought to the domain which was 
previously dominated by phenomenological approaches. 
The obtained results may have an important consequence 
for astrophysics, namely for the approaches to r-process 
nucleosynthesis: those involving Brink-Axel hypothesis 
may need to be revised. 
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